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ABSTRACT 


During the past two decades, naturally occurring stable 
isotopes of oxygen and hydrogen have been used in different labora- 
tories to study the precipitation formation process. Although this 
technique has been known for several decades relatively little progress 
has resulted due to field, laboratory and theoretical difficulties. 
In the case of hail studies, field sampling of precipitation and 
water vapour is hazardous and difficult to control. Laboratory 
analyses require specialized sample preparation techniques and 
sophisticated mass spectrometers. 

A program of sampling vapour and precipitation jin Alberta 
was undertaken during the springs and summers of 1974, 1975 and 1976 
in order to understand and model precipitation processes, typical 
of this region, with stable isotopes. A simple convective cloud 
model has been considered to interpret the observations. The 
temperatures of formation of hail and graupel were calculated 
employing this model. Results show that small hail and graupel 
appear to collect most of their mass in the warm regions of the 
cloud. 

Isotopic values obtained in the present study have been 
compared with established relations between 6D, 6180 and mean tempera- 
tures. Satisfactory agreements are obtained in most cases. Suggest- 


ions for further work have been made. 
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CHAPTER | 


INTRODUCTION 


1.1 General 


Most of our usable water is released in the form of precipita- 
tion from the atmosphere. For over a century scientists have attempted 
to understand precipitation and other associated phenomena. It is 
only during the last three decades that the broad subject of cloud 
physics has become better understood mainly because of technical de- 
velopments in the probing instrumentation. A detailed review of this 
subject is beyond the scope of this study. However, excellent books 
on cloud physics by Byers (1965), Fletcher (1969), Mason (1971, 1975), 


and Rogers (1976) are available. 


1.2 Isotopic Content of Water and Its Designation 


Atoms of many elements are found in nature having different 
atomic weights. Those atoms which have the same atomic number but 
different atomic weight are called isotopes. The chemical properties 
of isotopes of a given element are almost identical because the 


properties depend mainly upon the orbital electrons. The isotopy of 
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elements influences physical processes like diffusion, absorption, 
evaporation and condensation, the last two being of special interest 
to the present study. 

A water molecule consists of two hydrogen atoms and one 
oxygen atom. Hydrogen has two stable isotopes +H (Hydrogen) and 2H 
(Deuterium or D), whereas oxygen has three (160, !70 and 180). In the 
case of water, HD!60 and H,?80 are the most abundant after He 20, An- 
other isotope of hydrogen, 3H (Tritium or T) which is unstable is 
worthy of mention here. In the environment tritium comes mainly from 
thermonuclear explosions. Production of T by cosmic radiation is of 
minor significance. The exchange reactions between T and light water 
molecules give rise to tritiated water (HT0). These have been used 
as tracers to study the natural water systems (International Atomic 
Energy Agency, 1969, 1970, 1971). 

The abundance of ihe in natural waters is 2000 parts per 
million. The absolute concentrations of isotopic species are diffi- 
cult to measure accurately. However, the difference in the 189 content 
of two samples can be determined with high precision using mass 
spectrometric techniques. The measurements are reported giving delta 
values (6) defined by the expression 
Jive 02 a (parts: ber thousand)e (1) 


6(°/g9) = |(R 


een aoe renders 


where R and R are the 180/160 or D/H atomic ratios in 
sample standard 


the sample and standard, respectively. The reference standard employ- 


ed is SMOW (Standard Mean Ocean Water) as proposed by Craig (1961). 
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When a sample 6 value has been measured in terms of local or working 
standard, its value on the SMOW scale can be obtained from the 
following relationship, provided that 6 is known relative 


standard 
to SMOW. 


sample sample standard rs sample Standard 


6 6 6 6 x 10 3 97/99 (2) 


§ = 
SMOW standard SMOW standard SMOW 


1.3 ltsotopes in the Hydrological Cycle 


Isotopic fractionation of water, i.e., changes in the frac- 
tion of Hee or HD!®0 in nature can be attributed to a number of 
physical, chemical and biological processes. In the atmosphere, the 
following processes lead to changes in the abundance of isotopes 


(fractionation): 


1. Evaporation of vapour from water body surfaces and 
land and evapotranspiration from vegetation. 

2. Condensation of vapour during cloud formation. 

3. Isotopic exchange between vapour and its environment. 


4, Evaporation from falling precipitation. 


Comprehensive discussions of this topic have been given by Dansgaard 
(1964), Craig and Gordon (1965) and Dincer (1968). The main results 


of the previous studies can be summarized as follows: 


1. The isotopic composition of ocean water is nearly constant 
with concentrations of 997680:2000:320 ppm (parts per million) 


for H5+°0, H»180 and HD!®0, respectively. For instance 
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SMOW 
(Craig, 1961). 

The isotopic species are fractionated during evaporation and 
condensation processes because of differences in transport 
rates to and from water body surfaces, land and vegetation. 
When liquid and vapour are in equilibrium the ratio a of 

the atomic ratio of H» 180 to H, 160 in liquid (R ) to that in 
the vapour (R) is dependent on temperature (T) only. This 
ratio of the abundance of the heavy isotope in the liquid 
relative to that in the vapour is proportional to the ratio 
of the saturation vapour pressure of the common isotope 
(e169 


with Raoult's law. This ratio is called the fractionation 


) to that of the heavy isotope ( ) in accordance 


“H.180 


factor (a) and tt can be written as: 
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The fractionation factor decreases with increasing tempera- 


CUMS - 


A linear correlation between the mean annual surface tempera- 
ture T(°C) and 6!80 values was derived by Dansgaard (1964) 
from available data for high latitude non-continental atmos- 


pheric precipitation. This has the form 
§189 = 0.69 T- 13.6 (9/o9) . (3) 


It was also concluded that heavy isotope content of precipita- 


tion showed a minimum in winter and a maximum in summer. 
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>. On the basis of 400 samples of water from streams, lakes and 


Precipitation, Craig (1961) established the following linear 


relation between 6D and $180: 
éD a 1896450 + 10 (4) 


The value of slope and the constant depends upon the air 
humidity, kinetic processes during evaporation, and isotope 
exchange processes. The "'universal'' nature of this approx- 
imate equation has been confirmed by more recent studies 
(Rychkov and Vetshteyn, 1976; Hage, Gray and Linton, 1975; 


Miyake, Matsubaya and Nishihara, 1968). 


ee eHanl Studies 


The first application of isotopic technique to obtain in- 
formation on the growth environment of a hailstone was demonstrated 
by Facy, Merlivat, Nief and Roth (1963). These authors determined 
the D/H ratio of the successive layers of a hailstone and used these 
values, along with the knowledge of the vertical structure of the 
atmosphere during the shower, to trace its life history. They assumed 
that the isotopic content of each layer of the hailstone represented 
that of the supercooled cloud droplets which were accreted during 
its formation. During condensation, isotopic fractionation leads to 
enrichment of the liquid phase in heavy isotopic molecules, HDO and 
H,+80. As the condensation proceeds and the air parcel rises the 
Soe vapour becomes progressively depleted in the heavy isotopic 


species. Thus the isotopic analysis provides a measure of temperature 
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and consequently the height at which that layer was formed. 

Since then a number of articles on the isotopic analysis of 
hailstones from different regions of the globe have been published 
(Majzoube, Nief and Roth, 1968; Baily, Hulston, Macklin and Stewart, 
1969; Macklin, Merlivat and Stevenson, 19702 "Knight, chhalt, Roper 
and Knight, 1975; Jouzel, Merlivat and Roth, 1975; Macklin, Knight, 
Moore, Knight, Pollock, Carras and Thwaites, 1977; Ross, Schooling 
and Vogel, 1977). None of these studies includes measurement of iso- 


topic ratio of the vapour feeding the cloud. 


1.5 Modelling of Precipitation 


It is well known that a cloud is formed by adiabatic cooling 
of a saturated air parcel during ascent. The major process in the 
formation of rain, hail or snow involves the capture of cloud drop- 
lets. The isotope content in the precipitation depends upon the iso- 
tope content in the cloud droplets. When this capture involves the 
instantaneous freezing of cloud droplets, such as during graupel 
(snow pellet) formation, the 6 value of the cloud droplets is ''frozen'! 
into the particle for future field sampling. Liquid precipitation on 
the other hand tends to change its isotopic value during subsequent 
growth in the cloud and during evaporation below the cloud base. Two 
types of models for the cooling processes, namely closed and open models, 
under equilibrium and non-equilibrium conditions, have been described 


in the literature to explain the isotopic content of cloud and vapour. 
In a closed model, cloud droplets and vapour are not separ- 


ated during cooling, i.e., the cloud droplets are not precipitated 
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Out. In the case of an open model, some of the cloud droplets are 
precipitated out of the system. This causes the remaining vapour to 
become more quickly depleted in heavier isotopes than in the closed 
model. 

Non-equilibrium models attempt to model processes in which 
the delta of vapour relative to the delta of the condensed water is 
not determined by the fractionation factor during evaporation or con- 
densation. For instance, in Rayleigh distillation 6 is determined 


vapour 


by 6 and the fractionation factor: therefore, 6 is in- 


liquid vapour 
dependent of the vapour which had previously evaporated. 

Miyake, Matsubaya and Nishihara (1968) have studied the 
time variations of the isotopic content of precipitation of both warm 
frontal (continuous) and cold frontal precipitation (showers). With 
warm fronts the heavier isotope content in precipitation increases 
with time at a fixed point on the ground. In showers, however, no 
simple rule in isotopic variations with time is observed. This is 
because the delta value of precipitation is lower if it forms high 
in a cloud and higher with greater evaporation below the cloud base. 

The latest attempt in modelling hailstone growth using iso- 
topes is that of Jouzel, Merlivat and Roth (1975). They have analyzed 
and discussed variations in the isotopic content of hailstones on 
the basis of the model proposed by Chisholm (1970). According to 
this model, cloud is divided into two regions; an updraft core un- 
affected by entrainment phenomena and the boundary zone influenced 


by these phenomena. The relevance of the above models will be dis- 


cussed further in a later chapter. 
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1.6 Objectives of the Study 


This study was initiated with the following four objectives: 


|.) Tov obtain Ho 80/H, 160 and HDO/H,0 ratios in vapour for con- 
tinental air masses and relate this data to precipitation 


data on an hour to hour basis. 


2. To interpret and explain isotope variations in vapour and 
precipitation in terms of simple cloud models of showers 


and thundershowers. 


3. To use a simple cloud model to estimate the temperature at 
which snow pellets and small hailstones form since these 


particles are potential embryos for large hailstones. 


4, To compare oxygen-18/oxygen-16 variations in various forms 
of precipitation from Alberta with published studies from 


other regions of the globe. 
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CHAPTER 2 


EXPERIMENTAL PROCEDURES 


2.1 Procedure for Collecting Water Vapour 


Chapter 1 reviewed previous research on the isotopic con- 
tent of precipitation from convective clouds. This study was designed 
to monitor the isotopic content of the water vapour feeding such clouds 
as well as the precipitation which is formed. To perform an adequate 
number of experiments it was frequently necessary to collect water 
vapour samples from a travelling van (enclosed light truck). Porta- 
bility dictated the design of the vapour condenser. A forty kilogram 
block of dry ice was found to last for one week when stored in the 
van's freezer. This made extended trips possible during which vapour 
could be collected using dry ice as a coolant and to preserve hail- 
stones. 

Vapour was collected using standard condensers or Horibe 
traps made locally from glass tubing. Outer glass tubing (1.2 cm o.d.) 
was coiled around a center glass tube (5 cm o.d.) which had knobby 
restrictions to air flow to give the condenser an approximate total 
diameter and height of 8 and 30 cm, respectively. Air was sucked first 


into the broader center tube and then around the outer coil by connect- 
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ing the condenser with polyethylene tubing to a twelve volt air mattress 
pump. During suction, the condenser was submerged in a dry ice- 
ethanol slush bath. The ideal temperature of such a slush bath is 
-70°C, but it was found that -55°C to -60°C was easier to maintain in 
the field. Slush bath temperatures were maintained and noted during 
field experiments. The pumping speed of the twelve volt air mattress 
pump was lowered by reducing the cross-terminal voltage to about four 
volts. 
The equipment was mounted in the van during field experiments. 
The air intake was at the front of the van. It took between one and 
two hours to collect two to three milliliters of water-equivalent ice 
in the condenser. This amount of water was required for the oxygen 
isotope analysis. About one-half to one-quarter of a cubic meter of 
air passed through the condenser in the one to two hour sampling time. 
During the sampling time the ice would occasionally build up and reduce 
air flow at the point where the air left the large diameter center 
tube and entered the outer coil. At that time the pump was turned 
off, the condenser removed from the slush bath, and the bottom of the 
condenser was submerged in a tray of water to melt the ice blockage. 
The condensing experiment was then resumed. When enough ice was 
collected the condenser was removed from the slush bath. The air 
entrance and exits were plugged quite tightly with foam rubber, and 
the ice was allowed to melt and collect at the bottom of the condenser 
for a few hours so that it could be poured out into a bottle. 
Experiments were performed in the laboratory to determine 
the influence of the rate of air flow and slush bath temperature on the 


isotope content of the sample and on the fraction of the water vapour 
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condensed. The results of these tests were favourable with delta 
values generally within 0.5 parts per thousand and about 95 percent 

of the vapour was collected under typical conditions. These tests 
were mainly done in the third summer of sampling (1976) after previous 
experiments indicated that added confidence in the vapour collecting 
technique should be developed. 

By the third summer of experiments (1976), it became evident 
that the isotopic content of near ground water vapour varied consider- 
ably during typical summer days due to varying ground wetness and day- 
time atmospheric mixing. To get instantaneous vapour samples, some 
large plastic bags of volume approaching one-half of a cubic meter 
were made up. Using these bags to store air samples, it was possible 
to reduce the vapour sampling time from one or two hours to a couple 
of seconds. In We way breezes were monitored as frequently as seven 
times as a thundershower approached. The bags were stored in the van 
and the vapour condensed by pumping the air out at a later time under 


condenser conditions similar to earlier ''non-bag'' experiments. 


2.2 Procedure for Collecting Precipitation 


During the spring of 1974 experiments were designed to 
collect graupel (snow pellets) from convective clouds on the same 
dates as water vapour samples were collected. Commercial size plastic 
garbage bags were held open for a minute or two to collect the snow 
pellets or rain which was falling. Sometimes mixed rain and smal] 
hail were collected. A small bottle was used to quickly dip the 


precipitation from the bottom of the sampling bag before significant 
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evaporation could occur. Occasionally freshly fallen graupel was 
picked off the ground. Several samples were frequently taken from a 
single shower. A field trip to the mountains that spring resulted 

in snow also being sampled by the bag technique and by removal from 
the snow-pack. Most precipitation samples were about ten milliliters 
in volume. 

With the coming of summer, snow pellets seldom reach the 
ground on the plains before melting or at least partly melting. Two 
techniques were devised to separate rain from the unmelted ice pellets. 
The ice pellets were to be analyzed for isotope content for comparison 
to the vapour's isotope content so that one could estimate the cloud 
temperature at which they formed. The first technique for separating 
rain and hail used a box with sponge rubber in the bottom. When ex- 
posed to precipitation the rain water was soaked up by the sponge and 
the ice particles stayed on top. After bottling the ice particles 
the rain water was wrung out into another bottle. The second method 
of separating hail from rain had the precipitation coming through a 
funnel in the roof of the van after which the hail was deflected by 
a screen and the rain fell through the screen. 

Hail was frequently picked off the ground and stored in the 
van's freezer for later analysis and bottling. Frozen hail samples 
were also available from the Alberta Hail Project (ALHAP) sampling 
program and from individuals who had stored hailstones in their 
freezer. 

Sequential precipitation samples were frequently taken to 


get minute by minute changes during the showers. Bag samples were 
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generally not timed accurately enough to get precipitation rates but 
standard ALHAP beakers were often used to collect timed precipitation 


samples so that precipitation rates could be calculated. 


2.3 Conduct of Field Experiments 


During the first spring (1974) two hour vapour samples were 
collected prior to the onset of graupel showers. These showers were 
sought out using a car and the vapour samples were taken at the 
University of Alberta. Frequent weather notes were made during these 
experiments. The only field measurements besides precipitation were 
wet and dry bulb temperatures. These experiments were conducted to 
try to measure the temperature of graupel formation by using the iso- 
topic content of the vapour and the graupel particles. 

Late in May of 1974 a trip to the mountains was made to test 
the newly-equipped van and to develop sampling techniques for al] 
precipitation types. These experiments are discussed in the section 
on 1974 results. 

During the summer of 1974 the weather radar situated at 
ALHAP's research site near Red Deer was used along with a two-way 
radio in the van to guide the experiments to sites of probable hail- 
fall. Water vapour samples were taken during the morning and again 
on the way to shower activity. Vapour sampling was interrupted during 
showers since it would not be representative of air feeding the clouds. 
Frequent notes on paper and on tape recordings were taken during these 
"hail chases'! as wel] as during the two hours of vapour sampling. All 


sampling was done by experienced meteorology students. 
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During the summer of 1975 emphasis was put on water vapour 
sampling and seeking out heavy rain and heavy mixed rain and hail 
showers because our 1974 experience indicated that the influence of 
evaporation on the isotopic content of rain water could be eliminated 
during very heavy precipitation. It was unfortunate that even using 
radar guidance, little success was achieved that summer in locating 
very heavy precipitation events. These experiments were based on the 
theory that when heavy rain and hail fell simultaneously, most of 
the rain was from melted hailstones. In this case the isotope content 
of the rain, small hail and vapour would tell whether the unmelted 
haitlstone centers were formed at colder or warmer levels in the cloud 
than the rain which would come from melted outer layers of hailstones 
or melted small hailstones. 

During 1975 the two hour samples of water vapour were showing 
considerable diurnal variations. This indicated a need for instantaneous 
vapour samples and a study of the variation of the isotope value of 
vapour during the approach of thundershowers. 

The spring of 1976 was devoted to tests of the vapour sampling 
techniques under laboratory conditions. 

The summer of 1976 was devoted to frequent sampling of water 
vapour through the bag sample technique as showers approached and to 
the seeking out of heavy shower events using the 1975 techniques. 

The results of these three years of experiments are covered 
on a day-by-day basis and summarized in the third chapter of this 
study. It should be stressed again that extensive field notes were 


necessary to make sense out of these case studies. The objective of 
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measuring the isotopic hydrology of individual thundershowers was 


always the prime objective of the observers. 


2.4 Meteorological Data 


Surface wet and dry bulb temperatures were measured in the 
field using a sling psychrometer. The meteorological data needed to 
calculate the isotope model in the next chapter are given in Appendix 
A. Some of the information contained in this Appendix was obtained 
either from the Atmospheric Environment Service of Environment Canada 
or from the Alberta Hail Project. Some important weather observations 


taken from field study notebooks are also given in Appendix A. 


2.5 Oxygen Isotope Analyses Procedure 


Water poses serious problems when introduced into the mass 
spectrometer ionization chamber for isotopic analyses. However, 
189/169 ratios of water samples can be obtained indirectly by analyzing 
carbon dioxide equilibrated with a given liquid sample. The following 


exchange reaction occurs: 
H.180 (liquid) + C160, (gas) ZH21°O (liquid) + C60'80 (gas). (5) 


For precise determinations, the amount of C0, used must be considered, 
particularly if its isotopic composition differs significantly from 

the sample. Linton (1972) has discussed the consequences of ignoring 
this correction under experimental conditions similar to those employed 


in the present study. For less precise measurements, it setnot 
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necessary to know the fractionation factor a, defined as 


(189/169) 
CO> 


18 16 
(2°97 0) 1450 


The fractionation factor is temperature dependent and decreases with 
increasing temperature. If water is the dominant component of the 


system, then 


6120 = 618 
system Sneteen (7) 
Isotopic separation factor of C0.-water system has recently been dis- 


cussed by Horibe, Shegehore and Takakakuwa (1973). 


2.6 Equilibration of Water Samples with Carbon Dioxide 


The following procedure, developed by Epstein and Mayeda 


(1953), was used. It is also described by West (1972). 


1. About 2-3 ml of water were pipetted into each of six 25 ml 
equilibration flasks. A drop of HNO, was added to help in 
establishing equilibration. 

2. After closing, the flasks were attached to the vacuum line. 
Samples were frozen with liquid nitrogen and then the air 
was pumped out with a fore pump. 

3. A thaw and freeze procedure was carried out for about 5 
minutes in a dry ice-ethanol slush bath and at room tempera- 


ture, and the air removed from the samples with a mercury 
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diffusion pump. 

4. Commercially available C0, was introduced into the flasks 
with the pressure ranging from 41-170 mm Hg (5.3-22.7 kPa). 
Flasks were closed tightly and then removed from the vacuum 
line. 

5. Equilibration was carried out in a constant temperature bath 
maintained at 26°C. The flasks were left in the bath for 
30-70 hours duration. 

6. After equilibration, the flasks were again attached to the 
equilibration vacuum line and water samples frozen in a 
dry ice-ethanol slush bath. 

7. Using a liquid nitrogen bath, CO, gas from each sample flask 
was then transferred through two slush baths to breakseals 
or screw-type flasks. In the present study, mainly break- 
seals were used to avoid leakage. 

8. The breakseals were then sealed and labelled for mass 


Spectrometric analyses. 


2] isotopic: Data 


Experimental work was carried out in the Physics Department, 
University of Alberta, using the equipment and facilities of Dr. J. 
Gray. The mass spectrometric procedure described by Coleman and Gray 
(1972) was followed for measuring !80/!60 ratios. 

Measured 6180 values relative to the local standard and with 
respect to SMOW are given in Appendix B along with the volume of sample 


and the pressure of CO, used for equilibration. Four 6D values, 
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reported in Table B-1] were determined by Dr. H.R. Krouse of the Univers- 
ity of Calgary. Calculations of delta relative to SMOW were made using 
equation (2). No corrections have been applied for various pressures 


of CO. used in the equilibration of water samples. 


CHAPTER 3 


RESULTS AND DISCUSSION 


3.1 Closed System in Equilibrium Model 


An attempt is made to interpret the isotopic measurements 
employing this simple model. It is an adiabatic model of a convective 
cloud in which vapour is in isotopic equilibrium with the condensed 
water and none of the cloud water is precipitated out. Also, there is 
no mixing of the cloud parcel with the air surrounding the cloud. 

The following relationships can be written in terms of N, 
the number of moles of various isotopic species in the closed equi- 


librium system. 


ni8 + nis be yi8 (8) 
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nis 
Fraction condensed F = —~— (11) 
16 
Si 
Nis 
Initial isotopic ratio in Vapour Re a (12) 
ni6 
vo 


The subscripts c, v and vo refer to condensate, vapour and initial 
vapour, while the superscripts 18 and 16 denote Ho 20 and H,1®0 
molecules, respectively. Rearrangement of the above equations with 


the elimination of N's leads to the following relationship. 


R 
vo 


(l/a-- 1)F + 1 


Equation (13) gives the isotopic ratio of the condensed cloud water 
(R) in terms of the isotopic ratio of the water vapour (RO) entering 
the cloud base at the lifting condensation level (LCL). Both F and a 
are functions of cloud temperature and they can be calculated at 
various altitudes by means of a tephigram if the cloud base (LCL) 
temperature and pressure are known. 


In terms of delta values, equation (13) can be approximated 


by: 


= = 3 = hy 
Store (10° + bg) (a 1)F (14) 
where oF is the delta value of the condensing water and ann is the delta 
value of the vapour at cloud base. For the delta values obtained in 


. O ° 
our experiments, oye can be approximated at -25 /oo for use in the 
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right hand side of equation (14). Furthermore, if a is written as 


Xx 


then for a = 1.0135 at -20°C, x = 13.5. Values of a used in this 
Study earen 1.0135. 10123. 1-006, and 1.0101 for temperatures of --20, 
-10, 0, and +10°C, respectively (Dansgaard, 1964). Therefore, x 
varies from around 13.5 at -20°C to 10.1 at 10°C. These are normal 
cloud temperatures employed in the calculations. The 6 values of the 
condensing cloud water oe can, therefore, be approximated at various 


temperature levels in the cloud by: 


Or 
it 
Or 


a ur ORS ES dee (16) 


From the values of x at various temperatures one finds that 
x increases by about 1.25 units per 10°C of cooling starting with 
x= 10.) at 10°C. A tephigram gives F as a function of cloud tempera- 
ture. Calculating equation (16) for oe one finds that cs falls off 
aye UN Vase Elikel 0 hese for the first, second, and third ten degrees 
of cooling above the cloud base, provided that the cloud base tempera- 
ture is around 0°C. If the cloud base temperature is yar (warm for 
Alberta) 6 would fall off by about 2.7, 2.5, and 2.0 °/oo for the 
first, second and third 10, 0 of cooling above the cloud base. Snow 
pellets, and to some degree hailstones, form when cloud droplets with 
characteristic on freeze by riming on their surface. This indicates 
the model's ability to estimate the cloud temperature at which a snow 


pellet formed by finding its 6 value and She to #0.2 °/oo. As in- 
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dicated above, the accuracy of snow pellet formation temperatures cal- 
culated with this model is greater in a cloud with a cold base than a 
cloud with a warm base since - falls more quickly with temperature in 
a cold-based cloud. This assumes that a snow pellet (graupel particle) 
or a small hailstone forms by freezing accreted supercooled water 
droplets which do not change their 6 value when they freeze. It is 
also assumed that the snow pellets or small hailstones are able to 
maintain their 6. value as they melt on descent from the freezing level. 

Figure 1] summarizes schematically how the 6 value of vapour 
and cloud water varies with altitude in this model. Arrows indicate 
how the & values would vary if the cloud has a significant precipitation 
efficiency. In that case, the cloud water which is rich in H,!®80 near 
the cloud base precipitates out, leaving lower values at higher alti- 
tudes. This precipitating model was called the Open System in Equili- 
brium by Miyake et al. (1968), and was also used by Facy et al. (1963). 
It has been compared to clouds caused by slow, warm frontal accent but 
it probably has some similarity to a cumulonimbus shower with high 
precipitation efficiencies as well. This model indicates that the 6 
values of the cloud droplets at a given altitude in the cloud should be 
lower as the precipitation efficiency increases. 

In summary, Figure 1 shows that on at cloud base will be 
about 10-14 °/oo (equal to x) higher than ve and that Je will decrease 
by about 4, 3, 2 S/oon respectively, for each successive 10°C of cool- 
ing in a convective cloud with a base temperature of Ons b. will 
decrease somewhat less quickly with altitude in a cloud with a warmer 


base temperature. 
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An analogous study to that described above, for the fraction- 
ation of HD!®0 in a convective cloud could have been used to check the 
results of the Ho! 80 model. Unfortunately, at the time of these ex- 
periments, deuterium analysis facilities were not readily available 
on the University of Alberta campus. As mentioned in the previous 
chapter, four water samples were sent to the University of Calgary 
for analysis, and the results are discussed later. Both HD!®O and 
H2+80 models could be used to calculate the cloud temperature at which 
snow pellets or hailstones formed. Reviews by Dansgaard (1964) and 
Craig and Gordon (1965) indicate that the difference between the 
precipitation formation temperatures calculated by these two models 
would be an indication of the non-equilibrium processes at work in 
the cloud and during melting and evaporation below the cloud base. 
Such a study would be a useful continuation of the work reported 
here, however, it would probably not successfully isolate the con- 


tribution of the various non-equilibrium effects. 


3.2 Discussion of 1974 Results 


The delta values of 1974 samples as well as some meteoro- 
logical parameters are displayed in Figure 2. Some of these results 
have already been discussed (Charlton, Gray and Sandhu, 1975). The 
observations are discussed on an event basis emphasizing those situa- 


tions when samples of both vapour and precipitation were collected. 


MAY 15. Three samples of small hail (3 dots in the symbol 
for graupel or small hail), two samples of rain water and one of 


vapour were collected, all within 20 km of Edmonton. Arrows in Figure 
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2 show temperatures of the cloud water at the corresponding 6 values, 
calculated using equation (16), a tephigram, 80 and the temperature 
data given at the bottom of Figure 2. Comparing the 6 values with 
various temperatures indicates that the small hail (=5 mm in diameter) 


collected cloud droplets, on the average, at the -10°C level. 


MAY 17. Calculation of the cloud temperature at which the 
graupel grew was made difficult by large fluctuations in surface dew 
point (2.5 to -2.5°C). This variation in surface moisture, plus some 
surface temperature variation, gives a cloud base temperature between 
2°C and ale Surface observations from three airports near Edmonton 
indicate that a cloud base temperature of -5°C is most realistic. In 
that case the _graupel must have grown at temperatures between -20 and 
-30°C. Also shown in the May 17 data are the T vs. oe values for the 
possibility of warmer cloud temperatures indicating that the graupel 
grew on an average at -18°C. This exercise clearly shows the need to 
know cloud base temperatures accurately if the model is to be used 


to estimate growth temperatures of graupel or hail. 


MAY 23. This was a showery day in Banff National Park. A 
vapour sample was taken while ascending to Bow Pass (altitude about 
2.1 km) where falling snow and snow pack was collected. A vapour 
sample was also taken while descending after the precipitation had 
ended. The vapour sample with the high oF (=2273 °/o0) was taken in 
the evening during descent towards Calgary and remote from showers. 

Snow from the top three millimeters of the snow pack at 
Bow Pass was collected for comparison to the falling snow since it 


seemed to have fallen recently. There is a discrepancy in these 
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samples of 2 °%oo which is difficult to explain. The cloud base was near 
to the altitude of Bow Pass where the temperature and dew point were 
ESS Ge Using this information and equation (16) the growth tempera- 
tures of the snow were about -5°C for the falling snow and -15°C for 
the top 3 mm of the snow pack. However, there are two reasons why 
this model does not represent snowfall from what was probably primarily 
orographic snow. First, there is likely a significant precipitation 
efficiency in which case the 6 value of the precipitation would be 
lower than that calculated at a given temperature level (see Figure 1). 
The second problem lies with the fact that the growth process of the 
snow in this case appeared to be primarily deposition (the opposite 
of sublimation) since the crystals were only lightly rimed. The reader 
is referred to Miyake et al. (1968) for calculations involving an 
open system in the non-equilibrium and the use of fractionation factors 
applicable to sublimation-deposition processes. 

The snow pack compares well in delta value with the mountain 
creek water collected the previous day if some enrichment of creek 


water through evaporation before runoff is taken into account. 


JUNE 3. hy is higher than on previous days, probably because 
the vapour sample was taken in the afternoon when thundershowers were 
widespread. 6 for two rain samples which contained some small hail 
are similar, but considerable evaporation of the precipitation must 
have occurred, since their 6 indicates a mean growth temperature of 
only 22eC@eThe vapour sample may indicate Heo enrichment of the 


vapour through evaporation of precipitation. 
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JUNE 5. The high deita values indicate that rather strong 
evaporation was probably occurring. The showers were light. The 
variation of delta in these rain samples also indicates considerable 
variation in the amount of evaporation below cloud. Strong evaporation 
is also indicated by the fact that, without taking into account evapora- 
tion, the model of equation (16) would indicate that the rain formed 


near or below the cloud base. 


JULY 3. No showers occurred on July 3, but on the previous 
day, many hailstorms occurred between Calgary and Edmonton. A vapour 


sample taken at Red Deer showed high 6 (-22.5 °/oo) in the mid-morning. 


JULY 4. A vapour sample was taken before 1800h when thunder- 
showers began in the area. ov was considerably lower (5) Pa ° 7/00) 
than the day before indicating regional advection of air with lower Sy 
than the previous day. Violent thundershowers began north of Red Deer 
around 1800h and continued until after midnight as they moved east- 


ward. 


1900h. Three consecutive samples of mixed rain and hail 
were obtained in about five minutes as a shower approached the station- 
ary sampling vehicle. During this time the 6 value of the precipita- 
EUOnm cl atl ONleetl/= Dat Ohad °/oo, indicating decreasing enrichment due 
to below cloud evaporation as the below cloud relative humidity in- 
creased. This phenomenon has been studied by Hage et al. (1975), 


Miyake et al. (1968) and Craig et al. (1963). 


2000h. A sample of hail was collected at the side of the 
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road. It had fallen from a different thundershower than the 1900h 
samples, but it yielded a 6 value much like the previous mixed rain 
and hail samples. The hail was about 5 mm in diameter. The growth 
temperature of the hail as indicated by the model is about -3°C which 
is unrealistic. Apparently the Sy for the afternoon vapour sample was 
unrealistically low in the sense that it did not represent the vapour 


feeding the evening thundershowers. 


2100h. Two samples of rain from another shower continue to 


yield similar delta values. 


2400h. A very violent thundershower occurred during which 
four centimeters of rain fell in ten minutes. Small hail was separated 
from the rain and its 6 value was determined. In such violent showers 
little below cloud evaporation could occur as is indicated by the very 
low 6 value of the rain, the lowest of all our samples. In fact, the 
rain's 6 value was lower than that of the hailstone cores which were 
separated from the rain. The outer part of the hailstone is assumed 
to have melted. The conclusion is that the center of the hailstones 
were formed at a warmer level of the cloud than the water which made 
up the rain. This is one further example of hail embryos coming from 
the lower levels of a cloud. Jouzel et al. (1975) found similar re- 
sults in three Alberta hailstones. This was the experiment which led 


us to seek very heavy showers during the next two summers. 


The first vapour sample was taken over a period of two hours 
followed by the second sample which was taken over a period of two 


hours during which showers were developing in the area. Studies in 


ie ee petal rey; 183 At dwt — 
ett vy af ake —_ ates we 


ie aT we. Anke rue =. Sa allele 


vl 
— au LJ 


ad ma Vitw a9 bud 4> sulin Gis sil oa fweland se ry) aA te ! 
m : 
Pa y | peri eqit helt ban 


me 
Vabine’ BP bag ind Pe Bly 4) foe! verse aar +: wit b)lesleeet mn 


_ _ 
| bea Vor’. yuloape wte'g “Ts 


om ot og =< Tt ie > 104 Ve7T’ : a /9@, 1 , 
7 . 


, 


~~. 


el aii! Vie ene : ET FO ce Pee f ba 71 7 
- 
9G) @@ Dae nl Were, aed a 


~ a i ll a 
fe ats Ve J ise héaPiphee ry ww wmsili ' Pie ooeee A AeInG a 
@ ¥ 
Sone ise)! las! le it a rods bIAY alee vont onan 
. i ‘ 


‘\@e" OL a “Uaale, i J AAN Wj wi ide pet iw. a i> hin aha. iii wy 


pjew ani 7 baigateal : | } ' ennwe bie va 7 


; i CIpy: (} Lex wy: BMT eAVS aorg Kl Po, at Yea ote 
‘ah ie 
. mt - 
wets Veedh ll? owl an @ ere i’ tuie tr “lads nat gawd! iow wg Pay ai: 


7 
a 


ae 264 1 woayyi i ws 16 tee Le ea aiabior eg wenn’ Niet 


7 
pesl ip aNd Wh Wie an vil ‘ $feul jas off bas | ad & 


Chait appeinape As nett all *) pay Yo i>) Vetoew @ je enn) 
bio poison, High So via@pad: wm Bite el, chef lend 
A _ 
~ foie gn ye) ip-i9. Pieilet untae te ‘bevel “ 
We ak Pa 7 
7 _ ae alg ores a i at: peers ore way 
a | _ 7 . | 4 Apita Poteet alt vw ct 


é 7 va - 


30 


1975 and 1976, to follow, indicate that oy varies considerably with 
time. All precipitation samples plotted in Figure 2 were from differ- 
ent showers, so they are difficult to compare. The higher ps value 
indicates that the samples were formed near the freezing level of the 
cloud. Clearly, the oF values from the early afternoon and late 
afternoon are not high enough to represent the vapour feeding the 


evening showers if below cloud evaporation is not taken into account. 


JULY 9. Two vapour samples collected in the late morning 
and early afternoon of July 9 have very high values (-20 °/oo), 
apparently because the ground was wet. It had rained continuously the 
day before. It was a relatively hot day, surface temperature hee 
and showers developed in the evening. 6 values of hail samples are 
about 5 °/oo lower than the rain which contained some ice. The growth 


temperature of the hail as indicated by the model is about = 20 sce 


JULY 10. The delta value of the vapour sample (-24.9 ° 700) 
is a representative value for a sunny summer day in Alberta. The ground 


was dry in Edmonton and the dew point was a typical summer value of 


ie) 


Shin 


JULY 11. This was a day of continuous rain so that one 
expects the vapour to be in equilibrium with the rainwater. In this 
case the rain samples would have the same 6 value as the cloud drop- 
lets at the LCL calculated from the model. Of course, the LCL does 


not exist since this is not a day of convective weather. 
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3.3 Discussion of 1975 Results 


Figure 3 displays the results of the 1975 experiments. The 
same general trends as those found in 1974 observations are evident. 
Those dates when samples of vapour and precipitation are available 


are stressed in this discussion. 


JUNE 16. Calculated temperatures of hail formation are 
: fe) , : 
high (0 to -5°C). It is evident that the vapour sample taken in the 
early afternoon does not have a high enough delta value to represent 


the evening showers. 


JUNE 17. The vapour sample was taken during a rain shower 


so that the rain is nearly in equilibrium with the vapour. 


JUNE 18. Four consecutive vapour samples were taken between 
morning and afternoon. Morning fog dissipated early. Nearly overcast 
stratocumulus developed by noon. In the afternoon the cloud dissipated 
as humidity fell. The first vapour sample has an average delta value 
because it is in equilibrium with the wet ground. The second sample 
has a high delta value, probably because this was taken just as the 
last surface water was evaporating. The last surface water would be 
“very Richedln Hs °O% The third vapour sample indicates a normal delta 
value that has returned with continued development of stratocumulus 
cloud. The fourth sample was taken late in the afternoon as the dew 
point fell and the cloud dissipated. Its delta value is low because 


dry air of lower delta value is being advected into the area. 
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JUNE 20. A vapour sample was taken in the early afternoon 
when light showers were developing. The first precipitation sample 
consisted of mixed rain and hail from early in a shower. Its high 
delta value indicates that considerable enrichment of heavy isotope 
was occurring. Later during the shower, rain and hail were separated 
and they had similar delta values indicating a temperature of formation 
in the cloud of about 45 {e, Their similar delta values also indicate 
little below cloud enrichment due to evaporation if one assumes that 
the liquid portion is melted hail. One final rain sample near the 


end of the shower had an intermediate delta value. 


JUNE 21. The large difference between the delta values of 
the vapour is much like that which occurred on June 18. This also was 
a day which began with the ground wetted by the previous day's shower 
activities. Since each rain sample was from a different shower, the 
large variation in the delta values of the rain may be attributed to 
variations in the vapour's delta values as well as to differing degrees 


of below cloud evaporation. 


JULY 3. Unfortunately, the vapour sample for this date was 
lost in the laboratory. This was a day of heavy thundershowers and 
some walnut-sized hail was reported but none was sampled. In the 
early evening three samples of pea to grape-size hail were obtained 
from two different storms. Their relatively high delta values (-10.8, 
= 1205 sande) 565 °/o0) indicate that the vapour feeding the cloud must 
have been rich in !80 by about 5 °/oo. This condition appears to 
have occurred again on August 6 and August 28, two more hail days 


when no vapour samples were collected. 
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AUGUST 5. The delta value of the vapour (-18 °/oo) and of 
one sample of hail (+6.8 °/oo) are probably unrealistically high. A 
second sample of small hail had a near normal delta value (-17.6 °/oo). 
There were problems with the mass spectrometer during the analysis of 


this date's samples. 


3.4 Discussion of 1976 Results 


Observations for 1976 are given in Tables A-III and B-II|I 
(Appendix). Final results are shown in Figure 4. In general, trends 
similar to those of 1974 and 1975 are observed in the 1976 measure- 


ments. Some of the interesting situations are discussed below. 


JUNE 19. No vapour sample was collected on this day prior 
to the showers. However, the 6 values of the precipitation samples 
were measured and are shown in Figure 4. This figure indicates that 
the 6 value for vapour collected on June 18 and 21 is around -28 ° foo. 
This 6 value is used to calculate the cloud temperatures. Arrows shown 
in Figure 4 indicate that the rain and graupel were formed between the 


-8 and 0°C levels. 


JUNE 30. Four bagged vapour samples were collected, the 
last three on the top of a 15-storey building. The day had extensive 
instability but very little precipitation was reaching the ground. 
The very high delta value (-9.5 °/00) was recorded at a time when 
virga (evaporating precipitation) was very close to the collection 


point. Completely evaporated rain could give such a high value. 


a ae 
a a nee Res 


u 


rie th. Aan . 
Gatland a eihienairze:: pane shar 


POPee. Ane Pii-s soldat ol navie eu Oe! o!) al eegseeee 
seuere, iia fi » Sip l?-n w.1@ ere O8lgde® Teal ¢ cin 
oy Ode? are Sl ©" Yeerseain ereh bas ATET 1 eee an 
rai 60) BieedslidlG wis éeaki aici siweavedn 1 Genk a bevtaly 


SETH yo eh re Oh iva! a tee viet amily wt 
iolwaas imatamp hater oo | Tayanp dey i git | hedee | 
14AP wees ftd! reer pie Spy) pet ome) +» bm betes 
ava” Bs (phate 47 Ae ie ht ahs oer Did wile rrpyev “ap? Relpy ‘ 
heii 2 wore! wibebingets bwais. wis | Belye) oh oi Tt cola’ oat | 
oh) S-a4idy Peed) ovew fanny tom lee Sheds Afay ibal ¥ wuiel “a 
sieved. 9° bag \- : 


ae ehh Uremml 7 ie Mlk, | 

Pr stild Mm Yeragasd! @ ie coed ori: 7 Py He 
mabtwe Lah +P! prey ‘nll yale 
Ma SOT Rela: ely h iall 


35 


"soipnys 9/61 pue S/61 “4L61 
au} JO} ULe4 pue [ey *4noded 4O san ea 
eJ,ep obesene Jussoide4 soul, poysep ey, 


“9/61 Ul uoize3idioe1d pue anoded u} 09 1°H *t ounbly 
18H © jiey jews 40 jadnes5 V uley 9) MOUS %& joden /\ 
OL O€L OTL OO OT 06 OU O€l OSGI OVI OF LS 001 Se OE Ou (Qo) a 
O'%Z O0Z OSGI OI ODT OFT 0-02 O12 O€2 O02 OZI TI O'v2 pol SOL O9T (Gal Sa 
G2 OOF OG Ce EE 92 26 S Ol GIlS6 S¢ 19 v8 8S OS B89 (84/8)) & 
| [see] | [eed [seer | (eee) l | | 
We, SI vl el 8 | O€ Id 6l 8I at SI | ral 9/61 
| Alar anne | | HOMW 
— | - —- SE 
M c 
I Poe 
iN A Piss ae A | £ 
: /\ OAQ A /\ A Bee. i A - es 
= eee aw ee Pe ate Sa (SA el el I Sempra) ee ee ee ee oe eee eae ery 2 SS poe YE, a ee 
AA Ong AA A AN “if WOowd |x Ese: 
AY he See = — / MONS |x. IE oa 
ong N\ .01- Vv ¢ = 
v a= 
ae ON —- 0Z- o 
Pgs ee ae i 2 me 
“Ya t————9,¢+-"- —— 9.0 V—-—+———+————————— —-————— + — — — — — — — — — - - - a 
reais | — 5 Jel+=101 : SS 
ee | eon ee 0 ee eee St 
a g| Joll+=101 [ 
0) lz 
I ms 


36 


JULY 8. Seven vapour samples were collected. Measured 6 
values are typical of the vapour from central Alberta. hy varies by 
about 4 °/oo. The causes of this variation are not yet clear. Calcula- 
tions, using equation (16) and two sets of values for surface tempera- 
tures and ae were performed. These calculations suggest that the 
small hail formed at the Dec or the oaG level. If the suggested model 
is correct, there is clearly a need for the vapour feeding the cloud 


to be higher in delta value than that measured at the ground. 


JULY 14. Four samples of vapour taken in the laboratory over 
a period of three hours showed little variation in their delta values. 
This shows that the vapour sampling technique can give consistent 


results. 


3.5 General Discussion 


Calculated mean 6 values from the data given in Appendix B 
are -26.0, -23.8, -17.9 and -14.6 for vapour, snow, hail, and rain, 
respectively. The mean value of 5, is lower by about 11 °7oo than 
the mean 6 value for rain. O'Neill and Epstein (1966) found that Sy 
fell to values of between -20 and -24 °/oo during a California stormy 
period while calm weather conditions gave oe = =|4 ooo Other studies 
of the isotopic content of water vapour include those by Ehhalt (1974) 
who reported the vertical profiles of HTO, HDO and H20 in the tropo- 
sphere for some places in the United States. In most cases the heavy 
isotope content of the water vapour decreases as a function of the 
altitude. These observations imply that Ty measured at the surface 


would be lower on days of deep convection. The large variations and 
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mean values of CF found in this study are consistent with those found 
by Facy et al. (1963) at a continental location. In Alberta, summer 
convective precipitation comes mainly from locally evaporated water. 
Mean oF values and variations of oa for maritime locations are con- 
siderably lower as reported in the studies discussed above and in Craig 
and Gordon (1965). 

In general, the rain samples show higher 6 values than 
graupel or hail. This does not necessarily indicate that rain was 
formed at a lower level in the cloud since evaporation below cloud 
base would enrich the heavy isotope in precipitation. Studies of 
this phenomenon, which do not necessarily occur in isotopic equilibrium 
or by Rayleigh distillation, have been published by Hage et al. (1975), 
Miyake et al. (1968), Dansgaard (1964), Craig et al. (1963), Friedman 
et al. (1962) and Stewart (1975). 

Four samples of precipitation (Table B-1) were analyzed for 
D content in order to check the established relationship between 6D 
and 6180. Craig (1961) observed the following linear relationship 
between 6180 and 6D using four hundred samples of naturally occurring 


water. 
60a" hse SOME ON (2/60) (4) 


Dansgaard (1964) using least square analysis on monthly values from 
individual stations found slopes and intercepts that varied from 3.5 

to 8.0 and -30 to +30, respectively. Recently, for Edmonton precipita- 
tion, Hage, Gray and Linton (1975) have calculated the slope and inter- 
cept of equation (4) to be 7.8 and +0.4, respectively. Using mean 


values of 6180 = -19.2 °/oo from four 1974 samples, slope = 7.8 and 


* J - y 7 
> : 
& ahi ey > 
a 
yenenia no vel 7" 


, hs : : 
7 wv Hoke ry oa yO nye) Giebere wns rai ”" 
ore “ate a padi lns at iL i: oneal 


viet pi ee We bAmuleta pellore mld AP bab tediy oe econ alae 
, eee) | 


nw etnled. & delimit ié4e o vias. ene . btteong vel 
ray whe o> wha ' er? 4 fide : seat «ae , ave? . hed so! inte j 
D ) es 
| wa Cie mp TOSS yh het! 8 emo! at) «4 tea ‘aves! a Des wed 
| = 
7 a 
Vip ; »4 66 > 405 | nH! re ae py 1 y ; th re ’ ve aiven: ial r 
ayia rtau nc oh ai 45a Vi é 13) of Maelie&y (aaa Ata mn 
5 ' 


51 ey) : ‘ 7 7 wail vt) bye ; 1% User « ® s 1% ; iis ; ‘avis me 


te wf '¢ en q; a * ihe oa ; * er i wi od a 
“| HGabl saaaee by (seg? +h 4 
(a °ae 


$07? Vaavi ¢ ae hal { “ ae | t re ' iy ie De e int 
7 7 : f, a. is 
OF Sesdigh pileigalucler bwtalliorcn ¢% Vea @ aad i 
a] 7 % ’ _ 
GPU ETAT ONT OY Naleel peti lal ident CIPO) wdaad a bi 


» 
_ 


prt iS it WUTC Pa 1A celia! Wolrduri que paleo one orl, or 


ey | ual’) Jat atin ers ‘pe 


vin cipy ti ona oes) i nal eee, oie v : 
ae Era Sale - tha Raunt eae 7 
: er | 


38 


intercept = +0.4, the computed value of 5D is found to be -149 °/oo. 
This value is in excellent agreement with the mean value of -150 °/00 
found in this work. 

Using linear regression on the mean monthly delta values and 
mean monthly temperatures (T, ve for only the summer months at Edmonton, 


Linton (1972) obtained the following equation: 


Suiy S Wes We = Aide (may. (17) 


The slope and intercept of equation (17) has been discussed in detail 
by Hage, Gray and Linton (1975). There is no intention to discuss 
the climatology of isotopes in precipitation and repeat previous dis- 
cussions here. It is sufficient to mention that the isotopic data 
obtained in the present study are consistent with equation (17) and 
support Hage et al.'s (1975) argument that almost all warm season 
precipitation in central Alberta, shielded by high mountain ranges to 
the west, originates from water vapour produced by re-evaporation of 
precipitation over land (including transpiration). 

Recently Jouzel et al. (1975) have presented a detailed 
interpretation and discussion of the isotopic measurements made on 
three hailstones produced by a storm which occurred on August 7, 1971, 
in Alberta. The maximum and minimum 6D values as a function of the 
radius within the hailstones, taken from their Figures 5 and 6, are 
between -125 and -150 Choon Using these values and Hage et al.'s 
slope and intercept in equation (3), 6180 has been calculated, and it 
is found to vary between -15.4 and -19.2 °/oo. This is in good agree- 


ment with the majority of the 6 values observed for the hail and graupel 
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CHAPTER 4 
CONCLUSIONS 


4.1 Summary of the Data and the Model 


This study was designed to sample vapour, hail and graupel 
on showery days. The samples were to be analyzed for aD and the 
resulting data were to be used to understand and learn about precipita- 
tion formation menreseee Some difficulties in controlling field, 
laboratory, and theoretical studies arose, but they did not obscure 
some very interesting observations. During the three years of these 
studies, individual field experiments became more directed and con- 
trollable. For instance, studies of the variation of the vapour delta 
values on an hourly or even on shorter time scale became a necessary 
objective. Also, the study of rain during very heavy showers to 
eliminate the below cloud evaporation problem became an objective. As 
the time span it became apparent that graupel, i.e., snow pellets, 
constitute the only precipitation mechanism that comes very near to 
really ''locking in'' the isotopic content of the cloud droplets. 

The majority of the results suggest that small hail and 
graupel appear to grow in the warm regions of the cloud (T = 0 to -15°C). 


This temperature range, calculated using the simple model, is warmer 
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than what one would expect. Use of a more detailed isotopic cloud 
model which includes a finite precipitation rate (open models) would 
make this temperature range even warmer. The experiments showed that 
the delta value of water vapour is quite variable on an hour-to-hour 
basis during showery weather in Alberta. Therefore, delta values of 
vapour measured at the ground must be applied with caution to isotopic 
models of convective precipitation formation. Simultaneous measure- 
ments of deuterium and oxygen-18 would give better insight as well as 
a ''double check'' into the reliability of the numerical models used. 

During 1974, the observed delta values of vapour were in the 
range of -31.3 to -18.5 °/oo. The majority of these values are close to 
@2b °/oo. Reasonable explanations of these variations have been given 
in terms of the hourly and daily hydrological cycle. Mean delta value 
for the rain (-15.7 +/ 90) obtained in the present study was about 
4 °/o0 higher than the previous studies. Only 23 samples of rain were 
collected and analyzed in 1974. Such a small sample does not represent 
the average seasonal delta value of rain for Alberta because many 
samples were from light showers. Under such conditions, evaporation 
during rain falling would cause the delta values to rise. The 1975 
and 1976 data lead to similar observations. 

The present results support previous studies which suggest 
that almost all warm season precipitation in central Alberta originates 


from water vapour produced by re-evaporation of local water. 


4.2 Suggestions for Further Work 


The observed isotopic variations have been interpreted using 
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a simple cloud model and surface meteorological data only. It would 
be worthwhile to investigate and interpret the present data further 
using all available upper air observations and more detailed models. 

The present and previous investigations on isotopes in 
precipitation suggest that a more practical use of models would in- 
volve some simple and continuous method of sampling the water vapour 
in the air. It is suggested that new methods and procedures of sampling 
water vapour at the surface and aloft should be developed. 

It is easy to suggest that someone sample cloud droplets 
under a variety of weather conditions. Isotopic information on cloud 
droplets would be crucial in the selection of a model for practical 
applications. However, such observations have been suggested for years 
but no one has come up with a method to sample cloud droplets without 
having them partially evaporate. 

It is suggested that future studies should include simultan- 
eous measurement of deuterium and oxygen-18. 

Experiments on cloud seeding to suppress hail have been under- 
way in Canada and the United States for more than two decases. Al] 
the scientific evidence obtained so far is inconclusive in that it does 
not answer the question: ''Does cloud seeding work?'' The variations 
of H,180 and HDO in precipitation from seeded and nonseeded clouds has 
not been studied previously. It is felt that the study of seeding 
material (Agl) in conjunction with the isotopes in total precipitation, 
both rain and hail, from seeded and nonseeded clouds might yield clues 


to assess the effectiveness of hail suppression experiments. 
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APPENDIX A 


Meteorological Data for the Sampling Period 


Explanation of Sample Code 


YR  - Year during which sample was collected 
MO - Month of the year 
DA - Day of the month 


LDT - Local daylight time 


TYPE - Type of sample 


V - Vapour 

R Seenaat 

H - Hail 

G - Graupel 

RH = - Rain and hail mixed 
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APPENDIX B 


Isotopic Analysis of Collected Samples 


Explanation of Sample Code 


- Year during which sample was collected 
- Month of the year 


- Day of the month 


LDT - Local daylight time 


TYPE - Type of sample 


RH 


- Vapour 
= edi 
cetlae| 
- Graupel 


- Rain and hail mixed 
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